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Abstract. For a difference approximations of multidimensional diffusion, the truncated local limit 
theorem is proved. Under very mild conditions on the distribution of the difference terms, this the- 
orem provides that the transition probabilities of these approximations, after truncation of some 
asymptotically negligible terms, possess a densities that converge uniformly to the transition prob- 
ability density for the limiting diffusion and satisfy a uniform diffusion-type estimates. The proof 
is based on the new version of the Malliavin calculus for the product of finite family of measures, 
that may contain non-trivial singular components. An applications for uniform estimates for mix- 
ing and convergence rates for difference approximations to SDE's and for convergence of difference 
approximations for local times of multidimensional diffusions are given. 



Introduction 

Consider diffusion process X in R d defined by an SDE 

(0.1) X(t)=X(0)+ f a(X(s))ds+ [ b{X (s)) dW (s) , (eR + , 

Jo Jo 

and a sequence of processes X n , n > 1, with their values at the time moments ^, k E N defined by a 
difference relation 
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and, at all the other time moments, defined in a piece-wise linear way: 
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Here and below, W is a Wiener process valued in R d , {£ k } is a sequence of i.i.d. random vectors in 
M. d , that belong to the domain of attraction of the normal law, are centered and have the identity for 
covariance matrix. Under standard assumptions on the coefficients of the equations (|0.ip , ([0.2p (local 
Lipschitz condition and linear growth condition), the distributions of the processes X n in C(M + ,M d ) 
with the given initial value X n (0) = x converge weakly to the distribution of the process X with 
X(0) = x ([I]). Thus, it is natural to call the sequence the difference approximation for the 

diffusion X. 

Consider the transition probabilities for the processes X,X n : 
Px,t(dy) = P(X{t) e dy\X(0) = x), P£Jdy) = P(X n (t) e dy\X n (0) = x), t > 0,x G R d . 
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It is well known ([2]) that if the coefficients a, b are Holder continuous and bounded and the matrix b-b* 
is uniformly non-degenerate, then P Xjt (dy) = pt{x,y)dy. The function {pt(x,y),t G 6 M d } 

(the transition probability density for X) possesses the estimate 

(0.4) p t (x,y) < C(T)Hexp f- ^-^H ; t<2>,y G R d . 

The general question, that motivates the present paper, is whether any (more or less restrictive) 
conditions can be imposed on the coefficients a, b and the distribution of in order to provide that 
P x 1 t (dy) = p™(x,y)dy for n large enough, the densities p n possess an estimate analogous to (10. 4h 
and p n converge to p in an appropriate way. Such a question both is interesting by itself and has 
its origin in the numerous applications, such as nonpar ametric estimation problems in time series 
analysis and diffusion models (see the discussion in the Introduction to [3]), the uniform bounds for 
the mixing coefficients of the difference approximations to SDE's (see [1] and subsection 4.1 below), 
the difference approximation for local times of multidimensional diffusions (see [S] and subsection 4.2 
below) . 

In the current paper, we consider the question exposed above in a slightly modified setting. For 
the distributions P n , we prove the result that we call the truncated local limit theorem. Let us explain 
this term. We show that the kernel P n can be decomposed into the sum P n = Q n + R n in such a 
way that both Q n and R n are a non-negative kernels and 

(i) for Q n , its density q n exists, satisfies an analogue of (|0.4p and converges to p; 

(ii) for R n , its total mass can be estimated explicitly and converges to 0. 

The kernel Q n represents the main term of the distribution P n and satisfies the local limit theorem; 
the kernel R n represents the remainder term, and typically decreases rapidly (see statements (iii) and 
(iii') of Theorem 11.11 below) . Such kind of a representation appears to be powerful enough to provide 
non-trivial applications (see Section 2 below). On the other hand, the conditions that we impose on 
the distribution of in order to provide such a decomposition to exist are very mild; in a simplest 
cases, these conditions have "if and only if form (see Theorem 11.21 below). Our main tool in the 
current research is a certain modification of the Malliavin calculus. 

Let us make a brief overview of the bibliography in the field. Malliavin calculus have not been 
used widely for studying the properties of the distributions of the processes defined by the difference 
relations of the type (|0.2|) , (|0.3|) . The only paper in this direction available to the author is [4] , where 
rather restrictive conditions are imposed both on the coefficients (b should be constant) and the 
distribution of (it should possess the density from the class C d ). A powerful group of results is 
presented in the papers [3], [6], where a modification of the parametrix method in a difference set-up 
is developed. When applied to the problem formulated above, the results of [3], [6] allow one to 
prove that p n converge to p with the best possible rate 0(-^=). However, conditions imposed on 
the distribution of in [H] ,[■<?]) are somewhat more restrictive than those used in our approach. For 
instance, condition (A2) of [3] requires, in our settings, to possess the density of the class C 4 (R d ) 
(compare with the condition (B3) in Theorem 11.11 below) . 

The paper is organized in the following way. In Section 1, we formulate the main theorem of the 
paper together with its particular version, that is an intermediate between classic Gnedenko's and 
Prokhorov's local limit theorems. In the same section, we discuss briefly some possible improvements 
of the main result. In Section 2, two applications are given. In Section 3, the construction of the 
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partial Malliavin calculus, that is our main tool, is explained in details. In Section 4, the proofs of 
the main results are given. 

1. The main results 

1.1. Formulation. Let us introduce the notation. We write || • || for the Euclidean norm, not 
indicating explicitly the space this norm is written for. The adjoint matrix for the matrix A is 
denoted by A* . The classes of functions, that have k continuous derivatives, and functions, that are 
continuous and bounded together with their k derivatives, are denoted by C k and C^, correspondingly. 
The derivative (the gradient) is denoted by V, the partial derivative w.r.t. the variable x r is denoted 
by d r . The Lebesgue measure on R d is denoted by X d . For the measure /i on *B(R d ), [i ac denotes 
its absolutely continuous component w.r.t. X d . Any time the kernel P n is decomposed into a sum 
P n = Q n + R n , we mean that the kernels Q n , R n are non- negative; the same convention is used for 
decompositions of measures, also. In order to simplify notation we consider the processes defined 
by (|0.ip . (|0.2p and (10. 3D for t G [0,1] only. Of course, all the statements given below have their 
straightforward analogues on an arbitrary finite time interval [0, T\. 

Through all the paper, k is a fixed integer, k > 4. We denote e(/c) = K 2k+2 2 • 

Theorem 1.1. Let the following conditions hold true. 

(Bl) a G cl d+2)2 {R d ,R d ),b G C { b d+2)2 (R d ,R dxd ) and there exists (3 = (3(b) > such that 

(b(x)b*(x)v,v) R d > 7|M| 2 > x,v G R d . 

(B2„) E\\C k \\ K < +oo. 

(B3) There exist a G (0, 1) and bounded open set U C R d such that 

dfi ac a ^ H 

-aW ~ WU) U " ^ 
Then P n can be represented in the form P n = Q n + R n in such a way that 

(i) Ox tidy) = (it(x,y) dy and q n — > p,n — > +oo uniformly on the set [5, 1] x R d x R d for every 
<S €(0,1); 

(ii) there exist constants B, C, 7 > such that, for t G [0, 1], 



Q?ix,y)< 



Cri exp f-lM j y\\l\ \\x-y\\<tBn^ 
Ct 2 exp I —771 «+! ||x — y\ ) , \\x — y\\ > tBn ii + 1 



in addition, for every p > 1 there exists C p > such that, for t G [0, 1], x, y G 

qUx,y)<c p rl fi+ l| - , '" tfl 



t 

(iii) there exist constants D,p>0 such that R^ >t (R d ) < D[rr e ^ + e _pn '], x G R d , t G [0, 1]. 
If the condition (B2 K ) is replaced by the stronger condition 

(B2 cxp ) 3x > such that Sexp[x||^|| 2 ] < +00, 
then the following stronger analogues of (ii), (iii) hold true: 

(ii') there exist constants C, 7 > such that 

q?{x,y) <cHexp(-^-p^) , t€ [0,1], x,y£R d ; 
(iii') there exist constants D,p > such that R^ t iR d ) < De- pnt , x G R d , t G [0, 1]. 
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Let us formulate separately a particular version of Theorem 11.11 The most studied partial case of 
(|U.2|) is a = 0, b = J R d. In this case, X n (1) is just the normalized sum n~2 J2k=i an< ^ the limiting 
behavior of the distributions of such kind of a sums is given by the CLT. For the densities of the 
truncated distributions, the following criterium can be derived. We denote by P n the distribution of 

Theorem 1.2. The following statements are equivalent 

1. There exists no G N such that [P no ] ac is not equal to zero measure. 

2. There exists a representation of P n in the form P n = Q n + R n , such that 

(2i) Q n {dy) = q n (y) dy and sup yeRd 



- (27r)-2 e - 



-> 0, n — ► oo; 

(2ii) i/iere exist constants D,p > suc/t t/iai -Rn(M d ) < De~ pn ,n E N. 

The well known theorem by Prokhorov states that the given above statement 1 is equivalent to 
Li-convergence of the density of [-P n ] ac to the standard normal density (see [7], Theorem 4.4.1 for the 
case d = 1). There exist examples showing that, even while P\ <C X d , the density of P n may fail to 
converge to the standard normal density uniformly (see [7J, Ch. 4 §3 for the example by Kolmogorov 
and Gnedenko). The criterium of the uniform convergence is given by another well known theorem 
by Gnedenko: such a convergence holds if and only if there exists no £ N such that P no possesses a 
bounded density (see [7], Theorem 4.3.1 for the case d = 1). Theorem 11.21 shows the following curious 
feature: under condition of the Prokhorov's criterium, some exponentially negligible remainder term 
can be removed from the total distribution in such a way that, for the truncated distribution, the 
statement of the Gnedenko's theorem holds. This feature does not seem to be essentially new; one 
can provide it by using the Fourier transform technique, that is the standard tool in the proofs of 
the Prokhorov's and Gnedenko's theorems. We give a simple proof of Theorem 1 1 . 2 1 using the partial 
Malliavin calculus, developed in Section 3 below. This illustrates that the partial Malliavin calculus 
is a powerful tool that allow one to provide local limit theorems in a precise (in some cases, an "if 
and only if) form. 

1.2. Some possible improvements. In the present paper, in order to keep exposition reasonably 
short and transparent, we formulate the main results not in their widest possible generality. In this 
subsection, we discuss shortly what kind of improvements can be made in the context of our research. 

1. Difference relation (10.2[> is written w.r.t. uniform partitions {0 = t° < t\ < . . . }, = ^, k G 
Z + ,n £ N. Without a significant change of the proofs, one can prove analogues of Theorem 11.11 for 

the processes, defined by the difference relations of the type (|U.2[) with -% replaced by ^ • \Jt^ — t^ 1 
and partitions {i^} satisfying condition 

3c,C,d,D>0: lim inf -# Ik £ < t, (t* - tt' 1 ) £ -,— 
rw+oo n ' n n 



> td. 



(1.1) 



lim sup — jf 

n— *+oo Tl 



t n 5i (t n t n ) G 



c C 
n n 



<tD, te(o,i]. 



2. One can, without a significant change of the proofs, replace the sequence of i.i.d. random 
vectors in (|0.2p by a triangular array {£ n> fc, k < n} of independent random vectors, possibly 
not identically distributed, having zero mean and identity for the covariance matrix. Under such a 
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modification, condition (B2 K ) should be replaced by sup n fc £ , ||^ nj fc|| K < +00, and condition (B3) by 

jr.. lac 

(B3') 3a,r>0,x n GR d : ™ J > aI B(w) A d -a.s., 

here \x n ^ denotes the distribution of £ n &, B(x,r) denotes the open ball in R d with the centrum x 
and radius r. Also, the phase space for £ n £ may be equal R m with m> d (note that the case m < d 
is excluded by the condition (Bl)). 

3. Under an appropriate regularity conditions on a, b, Malliavin's representation, analogous to 
(|3.25|) . can be written for the derivatives of the truncated density of an arbitrary order with respect 
to both x and y. Thus, after some standard technical steps, one can obtain the following estimate, 
that generalize statement (ii') of Theorem ll.lt for a given k, I £ N, 

d k+l . d+k+i ( 7||y-x| |2 

jq t {x,y) < C k+ it 2 exp 



under (Bl), (B2 exp ),(B3) and a e cf+ fe+z+1 ) 2 {R d , R d ), b e C^ +k+l+1)2 {R d ,R dxd ). 

4. Theorem 13 . 1 1 provides the truncated limit theorem without essential restrictions on the structure 
of the functionals. For instance, one can apply this theorem in order to obtain a truncated local limit 
theorem for difference approximations of integral functionals, etc. 

5. Like the Malliavin calculus for (continuous time) diffusion processes, the partial Malliavin cal- 
cucus, developed in Section 3, can be applied when the diffusion matrix is not uniformly elliptic, 
but locally elliptic, only. However, the changes that should be done in the proof are significant; 
in particular, Theorem 13.11 is not powerful enough to cover this case. Thus we postpone the de- 
tailed investigation of this case (and more generally, the case of coefficients satisfying an analogue of 
Hormander condition) to some further research. 

6. In the present paper, we concentrate on the individual estimates for the densities q™ and do not 
deal with the convergence rate in the statement (i) of Theorem II. li The (seemingly) possible way to 
establish such a rate is to write the Malliavin's representation, analogous to (13.251) . for the limiting 
density p and then construct both the functionals X n (t) = f n and X(t) = f and the corresponding 
weights T* n and T*, involved into the Malliavin's representation, on the same probability space 
with a controlled L2-distance between (f n , T^ n ) and (/, T-^). Since the question about the estimates 
in the strong invariance principle for the pair (f n , T* n ) is far from being trivial, we postpone the 
detailed investigation of the rate of convergence in Theorem II .11 to some further research. We remark 
that the modification of the parametrix method, developed in [3], [6], provides, under more restrictive 
conditions on the distribution of {£fc}, the best possible converence rate 0(^=). 

2. Applications 
In this section, we formulate two applications of Theorem 11.11 



2.1. Mixing and convergence rates for difference approximations to SDE's. Under condi- 
tion (Bl) and some recurrence conditions, the process X is ergodic, i.e. possesses a unique invariant 
distribution ^ nv (see [8],|SJ). Moreover, an explicit estimates for the fl-mixing coefficients and for 
the rate of convergence of P x t = P(X(t) € -|-X"(0) = x) to Hi nv in total variation norm are also 
available. 

The processes X n , restricted to are a Markov chains. The following natural question takes its 
origins in a numerical applications: can the mentioned above estimates for the mixing and convergence 
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rate be made uniform over the class {X n ,n > l,X}? This question is studied in the recent paper 
[3J, see more discussion therein. 

In this subsection, we use the truncated local limit theorem (Theorem [TTTJ in order to establish the 
required uniform estimates. Denote, by || • \\ var , the total variation norm. Recall that the /3-mixing 
coefficient for X is defined by 

/3 x (t) = sup E»\\Pm,X(0)=x)-P(-\X(0)=x)\\ var ^ , (£R + , 



sGM+ 



t+s 



where 3* = o~(X(s),s G [a, b}), P(-|3q, X(0) = 0) denotes the conditional distribution of the process 
X with X(0) = x w.r.t. 9^, and 

B 1 f~)B 2 = 0, B 1 UB 2 = M. m 

B X ,B 2 g 5 

The /3-mixing coefficient {{3™(t),t G ^+} for X n is defined analogously. 

Theorem 2.1. Let conditions (Bl) and (B3) hold true. Suppose also that 
(B4) there exists Ro > and r > such that 

(a(x),x) R d < — r\\x\\, \\x\\ > Ro. 

(B5) there exists k > : E 1 exp[x||^||] < +oo. 

Then, for every process X,X n ,n > 1 there exists unique invariant distribution fj,i nv ,fj,f . More- 
over, there exist no G N, a function C : M. d — > IR + and a constant c > such that 

\\Pxt~ ^invWvar <C(x)e- Ct , t£-Z + ,n>n , \\P x ,t ~ Hinv\\var < C {x)e~ Ct , t G M + , 

n 

< C(x)e- Ct , t G -Z+, n > n , ^(t) < C(x) e - rf , t G R+. 
n 

Remarks. 1. The statement of Theorem 12.11 is analogous to the one of Theorem 1 [1]. The main 
improvement is that the conditions {D\) - (-D3) of Theorem 1 [3] are replaced by (seemingly, the 
mildest possible) condition (B3). In addition, Theorem l2.ll unlike Theorem 1 [3], admits non-constant 
diffusion coefficients b. 

2. The mixing and convergence rates established in Theorem 12.11 are called an exponential ones. 
If the recurrence condition (B4) is replaced by a weaker ones, then the subexponential or polynomial 
rates can be established (see Theorem 1 [1], cases 2 and 3). We do not give an explicit formulation 
here in order to shorten the exposition. 

2.2. Difference approximation for local times of multidimensional diffusions. Consider a 
W-measure [i on M. d , that is, by definition ([10], Chapter 8), a <r-finite measure satisfying the condition 

( 



(2.1) sup / w d (\\y - x\\)/j,(dy) < +00 with w d (r) 

xeK m ^||y-x||<l 



d= 1 

max(— lnr, 1), d=2 
r 2 ~ d , d>2 



Every such a measure generates a W-functional ([TO], Chapter 6) of a Wiener process W on M. d , 

rt dfi 



(2.2) V s * = V S '\W) = J %{W{r)) dr, 0<s<t. 
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For singular fj,, equality (12. 2p is a formal notation, that can be substantiated via an approxima- 
tive procedure with [i approximated by an absolutely continuous measures ([10]. Chapter 8). The 
functional p is naturally interpreted as the local time for the Wiener process, correspondent to the 
measure /i. 

Next, let the process X be defined by (jO.ip and satisfy (|0.4p . that means that the asymptotic 
behavior of its transition probability density as i ^ 0+ is similar to the one of the transition 
probability density for the Wiener process. Then the estimates, analogous to those given in [10] . 
Chapter 8 provide that the ^-functional of the process X 

(2.3) p s > t = p°>\X) = j^(X(r))dr, 0<s<t 

is well defined. We interpret this functional as the local time for the diffusion process X, correspondent 
to the measure /it. 

At last, let the sequence X n ,n £ N of difference approximations for X be defined by (I0.2p . (|0.3p , 
Consider a sequence of the functionals p n (X n ) of the processes X n of the form 

(2-4) = = £ fJxJ^Y < s < t. 

k:s<-<t 

— n 

In Theorem 12.21 below, we establish sufficient conditions for the joint distributions of (p n ,X n ) to 
converge weakly to the joint distribution of ((p, X). Thus, it is natural to say that the functionals p n 
defined by (|2.4|) provide the difference approximation for the local time <p defined by (|2.3|) . For the 
further discussion and references concerning this problem, we refer the reader to the recent paper [5]. 

We fix x E R d and suppose that X n (0) = X(0) = x. We denote T = {(s,t) : < s < t}. In 
order to shorten exposition, we suppose fi to be finite and to have a compact support. Together with 
the functionals (p n that are discontinuous w.r.t. variables s,t, we consider the "random broken line" 
processes 



C' t = ^"'" - (ns - j + l)ipn" + (nt - k + l)p n " se 



k-1 k 



n n I n n 

Theorem 2.2. Let conditions (Bl), (B2@), (B3) hold true. Suppose also that 
(B6) F n {x) > 0,x G M. d ,n > 1 and £ sup F n (x) -> 0, n -> oo; 

x£R d 

(B7) Measures ^L n (dx) = F n (x)\ d (dx) weakly converge to fi; 

(B8) limlimsup sup / w d (\\y - x\\)fi n (dy) ^ 0. 

oj.0 n-^+oo x <ER d J\\y-x\\<5 

Then (X n ,ip n (X n )) (X,tp(X)) in a sense of weak convergence in C(M + ,R rf ) x C(T,1R + ). 

Remarks. 1. The statement of Theorem 12.21 is analogous to the one of Theorem 2.1 [5]. The main 
improvement is that the condition A3) of Theorem 2.1 [5] is replaced by (seemingly, the mildest 
possible) condition (B3). 

2. Once Theorem 12.21 is proved, one can use the standard truncation procedure in order to replace 
the moment condition (B2@) by the Lyapunov type condition "3(5 > : i?||^|| 2+5 < +oo" (e.g. [TT] 
Section 5). 

3. For examples and a discussion on the relation between conditions (|2.ip and (B8), we refer the 
reader to [5]. 
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3. Partial Malliavin calculus on a space with a product measure 

For every given n £ N and t E [0, 1], the value X n (t) is a functional of £i, . . . , £„ and thus can be 
interpreted as a functional on the space (M rf ) n with the product measure fx n . However, under the 
conditions of Theorem 11.14 M may contain a singular component and therefore it may fail to have 
logarithmic derivative. Thus, in general, one can not write the integration-by-parts formula on the 
probability space ((M d ) n , (53(E d ))® n , /x n ). We overcome this difficulty by using the following trick. 
Under condition (B3), the measure (i can be decomposed into a sum 

(3.1) = a ■ hjj + (1 — a) ■ v, 

where nu is the uniform distribution on U . One can write (on an appropriate probability space) the 
representation for {^} corresponding to (|3. 1 P : 

(3.2) £ fc = E fe • rj k + (1 - e k ) ■ Cfc, 

where % ~ ttjj, ~ is, and the distribution x of e& is equal to Bernoulli distribution with x{l} = a. 
This representation allows one to consider the family £i, . . . , £ n (and, therefore, the process X n ) as a 
functional on the following probability space: 

(3.3) Q = (R d x {0,1} x R d ) n , J= («8(M d ) «> 2* ' 1 } ® <8QR d )) n , P = (ttj; x x x i/) n . 

Now, the measure -k\j has a logarithmic derivative w.r.t. a properly chosen vector field, and some kind 
of an integration-by-parts formula can be written on the probability space (O, 3~, P) (see subsection 
3.1 below). The Malliavin-type calculus, associated to this formula, is our main tool in the proof of 
Theorems 11.11 11.21 We call this calculus a partial one because the stochastic derivative, this calculus 
is based on, is defined w.r.t. a proper group of variables, while the other variables play the role of 
interfering terms. In this section, we give the main constructions of the partial Malliavin calculus, 
associated to the representation (13. 2p . 

3.1. Integration- by-parts formula. Derivative and divergence. Sobolev classes. Denote 

n = $7i x n 2 x ^3, 

n 1 = n 3 = (R d ) n , n 2 = {o,i} n . 

We write a point w G in the form to = (rj, e, £), where 

v = (m,---,Vn)e (K d ) n , e = (ei,...,e„) e{o,iy\ ( = (a, ...,&) e (M d ) n 

and ?7fc = (rjki, ■ ■ ■ ,md),(k = (Ckl, ■ ■ ■ , (kd) ■ in this notation, the random variables r] k > £ k,(k are 
defined just as the coordinate functionals: 

7? fc (w) = T] k , s k (uj) = e k , ( k (u) = Cfc, w = (rj,e,() G O. 

Denote by C the set of bounded measurable functions / on ft such that, for every (e, C) 6 ^2 x 
the function /(■,£, C) belongs to the class C°°(]R rf ) and 

esssup||[V^7((77, £ ,C))|| < +oo, j e N, 

where denotes the gradient w.r.t. variable r/. 

For / € C and A; = 1, . . . , n, r = 1, . . . , d, denote by dk r f the derivative of / w.r.t. the variable r]^- 
Also, denote by H the space R dxn considered as a (finite-dimensional) Hilbert space with the usual 
Euclid norm, and by {ek r , r = 1, . . . , d, k = 1, . . . ,n} the canonical basis in it: all coordinates of the 
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vector e kr are equal to zero except the coordinate with the index kr being equal to one. For a given 
functions tp : M. d — ► R and 6 n : (R d ) n — ► [0, 1], define the stochastic gradient D by the formula 

(3.4) [Df}(v,e,0=e n (C)Y,^(m)-[d k rf}(v,e,C)-e kr , / € 6. 

k,r 

This definition can be naturally extended to the functionals taking their values in a finite-dimensional 
Hilbert space Y (actually, in any separable Hilbert space, but we do not need such a generality in our 
further construction). Given an orthonormal basis {y{\ in Y, denote by Q Y the set of the functions 
of the type 

i 

and put for such a function 

Dy = J2[Dfi}®yi. 

i 

It is easy to see that the definitions of the class G Y and the derivative D do not depend on the choice 
of the basis {yi}. By the construction, D satisfies the chain rule: for any two spaces Y, Z and for any 
fi, ■ ■ ■ , fm E C y , F £ C°°(Y m , Z), m > 1, 

m 

(3.5) F(/!,...,/ m ) eC z and D [f(/ 1s . . . , f m )] = J2l d j F )(fh • • • , fm) ■ Df r 

3=1 

We denote D°f = f,D 1 f = Df. The higher derivatives D 3 \ j > 1 are defined iteratively: -D J = 
p .. . .. p (note that the first operator in this product acts on the elements of C Y while the last one 

3 

acts on the elements of 1) ® y ). 

Everywhere below, we suppose that U is an open ball B(z,r) (this obviously does not restrict 
generality). We define the function ip in (|3.4p by ip( x ) = r 2 — \\x—z\\ 2 . Due to this choice, if) S C QC (M. d ) 
and if) = on dU. These properties of imply the following integration by parts formula: 

[d r f)(x)^(x)dx = - [ [d r ^](x)f(x)dx, feC 1 (R d ),r = l,...,d. 
Ju 

As a corollary of this formula, we obtain the following statement. 

Proposition 3.1. For every h £ H and every / S C, the following integration-by-parts formula holds 
true: 

(3.6) E(Df, h) H = -e( p , h) H f, p = e n (() ■ Y^rMim) ■ e kr - 

k.r 

The formula (|3.6p allows one to introduce, in a standard way, the divergence operator corresponding 
to the derivative D. For g £ Q H ® Y ; pu t 

(3.7) 5(g) = ~y~] [(p, e kr )gkri + (F>g kr i, e kr ) H ■ yi, g kr i = (g, e kr ® Ui)h®y- 

k ,r,l 

By the choice of the function ip, 5(g) G Q Y as soon as j 6 The chain rule (]3.5p and the 

integration-by-parts formula (|3.6p imply that the operators D and 5 are mutually adjoint in a sense 
of the following duality formula: 

(3.8) E(Df, g) H = Ef5(g), f £ e Y , g £ Q Y > H . 
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Since, for every p > 1 and every Y, Q Y is dense in L p (Q,P, Y), the duality formula (|3.8p provides 
that, for any p > 1 and any Y, the operators D, 5 are closable as densely defined unbounded operators 

D : L P (Q, P, Y) -> L p (n, P,H®Y), 6: L P (Q, P,H®Y)^ L P (Q, P, Y). 

Definition 3.1. The Sobolev class WI^iY) (p > l,m£ Z+) is the completion of the class w.r.t. 
the norm 

i 
p 



p,m 



j=0 



< +0O. 



Since D is closable in L p sense, there exists the canonical embedding of W p m {Y) into L p {Q, P,Y). 

We denote W™(Y) = C\ mp W p n2 (Y). If Y = R then we denote the corresponding Sobolev spaces 
simply by W™. 

3.2. Algebraic relations for derivative and divergence. Moment estimates. Let us introduce 
some notation. We denote by C a constant such that its value can be calculated explicitly, but this 
calculation is omitted. The value of C may vary from line to line. If the value of the constant C 
depends on some parameters, say m,d, then we write C(m, d). The latter notation indicates that 
the value of the constant does not depend on other parameters (for instance, n). If, in a sequel, the 
constant C is referred to, then we endow it with the lower index like Ci,C2, etc. We use standard 
notation {5jj~,j, k £ N} for the Kroeneker's symbol. 

For an H <g> H <g> Y-valued element K, we denote by K* the element such that 

(K* ,h® g ®y)H®H®Y = (K,g ® h (g> v)h®h®y, h,g £ H,y eY. 

For an X® Y-valued element g\ and X®Z- valued element g2, we denote by (gi,g2)x the Y®Z- valued 
element 

(gi,92)x= ^2 {gi,x h ®y h ) X ®Y ■ {92, Xh ® z h )x®z ■ [yi 2 ® z h ], 
h,h,h 

here {xi},{yi}, {z{\ are orthonormal bases in X, Y and Z, correspondingly. We also denote for an 
Y (g) X-valued element g\ and Z (g) X- valued element gi 

{gi:92)x = ^2 {9l,Vh ® x h)Y®X ■ {92,Zl 3 ®X h )z<S>X ' [Vl 2 ® z h\' 

Although the same notation (•, -)x is used for two slightly different objects, it does not cause misun- 
derstanding further. 

Consider the £(.£f)-valued random element (i.e., random operator in H) B, defined by the relations 



(3.9) [Be kiri ,e k2r2 j ^ = - \D(p, e klTl ) H , e k2 r 2 ) H = ~ d n{C) ■ h ± k 2 • [d ri d r2 i/)]{r) kl ) ■ V(%), 
k\ t 2 = 1, • • • , n, ri 5 2 = 1, • • • , d. We define the action of B on H ® Y-valued element g by 

Bg = }Xg, e k , r ® Vi)h®y ■ [Be k ^ r ] ® yi . 

k,r,l 

Using representation (|3,7p . one can deduce the commutation relations for the operators D,5, 
analogous to those for the stochastic derivative and integral for the Wiener process (the proof is 
straightforward and omitted; for the Wiener case, see [12], §1.2). 
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Proposition 3.2. I. If f £ G, g £ Q H ® Y , then f-ge Q mY and 

5(f-g) = f-5(g)-(Df,g) H . 

II. If g G Q mY , then 

£>[<%)] =Bg + d([Dg)*). 

HI- If 91, 92 G C H , then 

{D[S(gi)],g 2 )H = (Bgi,g 2 )H + ,92)h) + ([Dgi}*,Dg2)H®H- 

The main result of this subsection is given by the following lemma. 
Lemma 3.1. Let m,l G N,g G Wf^+^IT). T/ien t/iere exisfc 5(5) G W^ OT and 

(3.10) ||%)||2m,Z < C(m,l,d,1p) • \\ghm,2m+l-l- 

Remarks. 1. On the Wiener space, the typical way to prove estimates of the type (|3,1U|) is to 
use Meyer's inequalities for the generator L = 5 ■ D of the Ornstein-Uhlenbeck semigroup (see, for 
instance, [12] §2.4). Moreover, on the Wiener space, (|3.10p can be made more precise: the similar 
inequality holds with 2m + 1 — 1 replaced by I + 1. In our settings, it is not clear whether the operator 
5 ■ D provides the analogues of Meyer's inequalities, since it does not have the specific structural 
properties of the Ornstein-Uhlenbeck generator (such as Mehler's formula, hypercontractivity of the 
associated semigroup, etc.). Thus we prove (|3.1Up straightforwardly by using an iterative integration- 
by-parts procedure. 

2. Throughout the exposition, the function ^ is fixed together with the set U = B(x, z). However, 
when the constant C depends on the values of ip or its derivatives, we indicate it explicitly in the 
notation for C. 

In order to prove Lemma 13.11 we need some auxiliary statements and notation. For g G Q Y and 
m G Z_|_, we define the random variable \g\ m by 



^ 1 1 &9\\ "h® i® Y 

Lemma 3.2. If g G Q H ® Y then Bg G Q H ® Y and, for every m G Z + , 

\Bg\m < C(m,d,ip) ■ \ dim- 
Proof. Write Bg in the coordinate form: 

Bg= ^ (#> e fcri ®Vi)h®Y ■ h^r* ' [ e fcr 2 ® Vli 

where b kr \ r 2 = [Be kr \,e kr 2)u (recall that (Be^i r i,e^ r j)ff = as soon as k 1 ^ k 2 ). Write the 
Leibnitz formula for the higher derivatives: 

x 



D m (Bg)= H E {D* @ {9,e kT ,® yi ) H ®Y,®e k ^ 

k,r x ,T 2 ,l Oe2{ 1 " "• m > ki,ri,...,h m ,r m ieO 



(3.11) x [D m - #0 b krl 



>kr l r 2 1 K/\) e kiTi 



H ®( m -#e) 

ige i=i 



ekin ® e kr 2 <g> yi 
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where #0 m denotes the number of elements in the set 0. We write 

5*6= ^2 ( K D * @ (9,e k r^yi)H()SY,<^)e kir 

k^r 1 ,r 2 ,l,ki,ri,...,k m ,r m iS0 



X 



X 



jym #©5 fcrl ^ 2) 



kr x r 2 > e ki 



#®(m-#e) 

igO 8=1 



and estimate ||»Se||.H-<s(m+i)(«y. The function -0 belongs to C°° and is bounded together with all its 
derivatives on U. Thus, one can deduce from the representation (13, 9p and formula (13, 4D that 

\\D M b <C(M,d,tp), MeN. 

In addition, due to (|3.9|) . 

(D-#^ fcrlr2 ,(g)e te )^ Me) =0 

i£0 

as soon as fcj 7^ for some i ^ ®. Using these facts, we deduce that 
fci6{l,...,n},ri6{l,...,d},i0e i^O 



2 < 



Thus 

\\Se\\ H 0( m +i)^ Y 

<C(m,d,ip) [D* (g,e kT i <g> Vi)h®Y,&) efcrt) #e < 

fc,r 1 ,r 2 ,«,fc 4 e{l,...,n},r l e{l,...,d},ie© 260 

< C(m,eJ,V>)M#e ^ c OM,V>)Mm- 

Taking the sum over 

9 £ 2< 1 '-» m > and using the Cauchy inequality, we obtain the required statement. 
The lemma is proved. 

Proposition 3.3. I. Let g x G Q X ® Y ,g 2 G Q x ® z , then (gi,g 2 )x G Q Y ® Z and 

\{gx,92)x\m < C(m)|5i| m |0 2 |m, > 0. 
II. Lei 5 G C y and ^4 G £(F, Z), then Ag G C z and \ Ag\ m < \\A\\ ■ \g\ m ,m > 0. 



The second statement is a straightforward corollary of the chain rule (|3.5p . The first one can be 
proved using the Leibnitz formula; the proof is totally analogous to the one of Lemma 13.21 and thus 
we omit the detailed exposition. 

Remark. Taking X = R, we obtain that, for gi G C Y ,g 2 G C z , g\ <g> g 2 6 e y ® z with |#i (8) # 2 |m < 

C(m)|flfi| m |5 2 |m> m>0. 

Using iteratively statement II of Proposition 13.21 we obtain that, for g G and m > 1, the 
derivative -D m [o~(<?)] can be expressed in the form 

£ m [%)] =^ m (5)+5(G m (<7)), 

where F m (g) G G H ® m ,G m (g) G g- f/ ® (m+1) axe defined via the iterative procedure 

F (g) = 0, G (g)=g, G i+1 (g) = [DGiig)]* , F i+1 (g) = DF t {g) + BG t (g), i > 0. 

The mapping if i— ► if* is an isometry in H ® H ®Y , thus statement II of Proposition 13.31 provides 
that 

(3.12) \G m {g)\j < \g\ m+ j, m,j > 0. 
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Using Lemma E21 we deduce that 

(3.13) \Fm{g)\j < C(m,j,d,ip) • \g\ m +j, m,j > 0. 

Thus, in order to prove inequality (|3.10p for g G Q H , it is sufficient to prove that 

(3.14) E\\5(g)\\ 2 Y m <C(m,d,^E\g\ 2 2 Z-i 

for any m > 1, g G Q H ® Y and arbitrary Hilbert space Y . In order to prove estimate (|3.14p we embed 
it into a larger family of estimates. Consider the following objects. 

1. Numbers ko, ■ ■ ■ ,k v £ Z + (i> < 2m) such that fco + ■ • • + k v = 2m. Denote Ij = [ko H h fcj-i + 

1, fc H h kj] n N,j = 1, . . . ,v, I = [1, ko] n N (if k = then I = 0)- 



2. Function a : {1, . . .,2m} -»■ {1,... ,m} such that ^cr" 1 ^}) = 2 and # Jj ncr _1 ({i}) 
every i = 1, . . . , m and j = 0, . . . , v. 

Lemma 3.3. Let f £ G Y ® ko ,gj £ Q H ® Ym ' . Then 

V 

(3.15) E Y, (50,02/^)^0 ■n <y ((»'®fia W V®*i) <C(w,d,V)-^ 
h,...,i m i&h i=i ieJj 

Remark. The left hand side of (|3.14p can be rewritten to the form 



< 1 for 



lffo|« |flj|t>-i 



(3.16) E Y ^(5,y/Jy)5((5,%)F)^(5,y/Jy)5((^y/Jy)...5((5,yz m )y^((5,^ m )y). 

If t> = 2m, /co = 0, fci = • • • = kim = 1, 9o = 1) gi — • • " = 52m = 3 £ gff®y ? then the left hand side of 
(|3,15p coincides with the expression written in (|3.16p . Thus Lemma 13.31 implies estimate ()3.14p . 

Proof of the lemma. We use induction by v. For v = 0, conditions imposed on a can be satisfied 
if m = 0, only (i.e., if go is a function valued in R). Thus, for v = 0, (|3.15p is trivial since go < 
l<?o| = | <?o I o - For v = 1, conditions imposed on <r imply that Jq = {1, . . . , m},I\ = {m + 1, . . . , 2m} 
and the function a, restricted to either Io or Ii, is bijective. Thus the left hand side of (|3.15p can be 
rewritten to the form 

m m 
!i,...,2 m i=l i=l 

where 7r is some permutation of {1, . . . , m}. Using duality formula (13. 8p . we rewrite this as 

EfDgo^g) , 

where the operator £ L{H ® y® m ) is defined by 

(3.17) A[h®y h ®---® yij = h® Vl<1) ® ■ ■ ■ ® y Ijr(m) . 

One can easily see that ^ is an isometry operator, and thus Proposition 13.31 provides that (|3.15p 
holds true for v = 1 with C(l, d, ip) = 1. 

Suppose that, for some V > 2, (I3.15P holds true for all i> < V — 1. Let us prove that (|3.15p holds 
for v = V, also. For every l\, . . . , l m , take 

V 

g = (fll,® W <rW )y®*u / = (0O»® W <rW )y« fc o • n i (fe'®W.(i))y»'i) 
ie/i iei"o 
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and apply duality formula (13.80 . Then the left hand side of (|3.15p transforms to the form 



{i,...,2 n 



»eio 



H 



i=2 ie/j 



+ S (50,® y« CTW Wo • (-D (*((fl r r,®W a(i) )y»*r) J ) (fll 5 ®l/i ff(0 )y»*i x 

r=2 h,...,l m iElo \ \ iel r / ie/i / H 



(3.18) 



x II 5 ((«j.®J/k(o)y®*3 
ie{2,...,v}\{r} ie/j 



Let us estimate every summand in (13.180 separately. The idea is that every such summand can be 
written as 



(3.19) 



■rpfe® 



with v = V — 1 or v = V — 2 and some new m, feo, ... , g~o,g v , and thus the inductive supposition 
can be applied. 

Consider the first summand. Denote by J the set of indices r £ {1, . . . , m} such that <7 _1 ({^}) C 
Jo U ii. In order to shorten notation, we suppose that J = {1, . . . ,#J} (this does not restrict 
generality since one can make an appropriate permutation of the set {1, ... , m} in order to provide 
such a property). Take permutations ttq : Iq — > Iq and tt\ : I\ — > Ji such that 

[<T0 7ro](i) = i, ie {1, ...,#J}, [cro7ri](i) = i- fe , £e {feo + l,...,feo + #</}. 

Then the first summand in (|3.18|) can be rewritten to the form 



E E 



fc fco+fel 

^ I D(A^ go,§§ yi^ o{i)) )Y^o,( A 7n9l, ® ^ ( . l(l) )y® fe i 
h,—,l#j \ *=1 i=fco+l 



if. 



V 



x II ®^w)y®' 
i=2 ie/ 3 - 

where the operators A,r and A 7Tl are defined analogously to (|3.17|) . Denote 



then 



9o — {E l A 7T0 g ,A 7Tl gi) H(S)Yei ( # j), 



k k +ki 

Y ( ^(Ar so,®yi^ o(i)) )y®*o,(Ari5i, ® 



!i,...,2 



1=1 



i=fc +l 



if. 



5o, 



yi 



CT(7T (i)) 



i=#J+l 



fcl 



""(""lM) 



yfe +fc 1 -2#J 
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Put m = m — #J, ko = ko + ki — 2#J, k\ = ■ ■ ■ , ky-i = ky and let {Iq, 1%, . . . , Iy-i} be the 
partition of {1, ... , 2m} corresponding to the family {ko, • • • , fey-i}. Put g\ = g 2 , ■ ■ ■ gv-i = 9v (§0 
is already defined). At last, define function a by 



cr(t) 



a(vr (i + #J)), i=l,...,k -#J 

ff(7n(i + 2#J)), i = k - # J + 1, . . . , k + ki - 2# J 

a(i + 2#J), i = ko + k\, . . . ,2m 



Under such a notation, the first summand in (|3. 18)) has exactly the form (|3.19p with v = V — 1, and 
the inductive supposition provides that this summand is dominated by the term 



C(V-l,d,ij>) ■ E 



v-i 

lffolv-1 H \9jW-2 



Since A 7TQ ,A 7Tl are isometric operators, we can apply Proposition 13.31 and obtain that 

\9o\v-l = \(DA no go,A 7Tl g 1 ) H(g)Ye)(# j ) \ V - 1 < C(V-l)\DA no g \ v ^ 1 \A ni g 1 \ V - 1 < C(V-l)\g \v\gi\v-i- 

For every j = 1, . . . , V — 1, |^j|y_2 = |g./+i|y-2 < \gj+i\v-i- Thus, under the inductive supposition, 
the first summand in (13.181) is dominated by the expression given in the right hand side of (13. 15f) . 

All the V — 1 summands in the second sum in (13.180 have the same form and can be estimated 
similarly; let us make such an estimation for r = 2. Using Proposition 13.31 we rewrite this summand 
to the form 



li,...,l r . 



iElo 



iElr 



ieh 



H 



i=3 



ieii 



+E (#0,® ^ CT(i) Wo<M (^((52,® y/ <T(i) )y®*r),(5i,®yi ff(i) )y(8fc 1 J X 

Zl,...,Z m ielo \\ iel r ieh / H/ 



J=3 



yz f 



+E Y (so,® y^Wo ^ffe,® 2/z CT(i) Wr) , (£>Si,®y^ 



il,...,i T , 



ielr 



ie/i 



(3.20) 



i=3 



Let us show that, after an appropriate rearrangement of the indices i, every summand in (I3.20p 
can be rewritten to the form (|3.19|) . Such a rearrangement can be organized in the way, totally 
analogous to the one used before while the first summand in (I3.18P was estimated. Therefore, in 
order to shorten exposition, we do not write here an explicit form for the permutations of the indices, 
used in such a rearrangement. The first summand in (|3.2(jp has the form (|3.19p with v = V — 2, 
9j = 9j+2,j = ,V- 2, 



9o = ( A Wl [Bg 2 ] , A n2 [#i <g> g 
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where J 1 is the set of such i E {1, . . . ,m} that cr _1 ({i}) C /0U/1UI2 (we do not write here an explicit 
expressions neither for the permutations 7Ti, 7T2 nor for ko, ■ ■ ■ , ky~2, &)■ Under inductive supposition, 
this summand is estimated by 

v V 
C(V - 2,d,ip)E\g \v-2 JJ \9j\v-3 < C(V, d, tp) E\Bg 2 \v-2 ■ \91\v-2 • \g0\v-2 JJ \9j\v-3 < 

3=3 j=3 

V V 
< C(V,d,tp)E\g 2 \v-2- \gi\v-2- \90\v-2 IJ \9j\v-3 < C(V,d,iJj)E\go\v bjlv-l, 

j=3 3=1 

here we used Proposition 13.31 and Lemma 13. 21 The second summand in (|3.20|) has the form (|3.19|) 
with v = V - 1, gj = g j+ i,j = 2, . . . , V - 1, go = go, 

gi = [DlA^^A^gi 



where J 2 is the set of such i £ {1, . . . , m} that cr — 1 ({i}) c 1% U /2- This summand again is estimated 
by 

v 

C(V - l,d,ip) E\g \v-i\gi\v-2 Y\_\9j\v-2 < 

3=3 

V V 

< C(V,d,ip) E\g \v-i\\Dg 2 \v-2 • \gi\v-2- |flo|v-2 JJ \gj\v-3 < C(V,d, ip) E\g \ v JJ \gj\v-i- 

j=3 3=1 

At last, the third summand in ()3.20p has the form (|3.19|) with v = V — 2, gj = gj+ 2 ,j = 1, ■ ■ ■ ,V — 2, 

9o = (A ni g , ([DA n2 g 2 ]* , DAir 3 g 2 )H®ii) 

and again is estimated by 

v 

C(V -2,d,i;)E\~g \v-2\Yl\gj\v-3 < 

3=3 

V V 
< C(V,d,1>)E\go\ V - 2 \ ■ \Dg 2 \ v _ 2 - \Dg x \ v _ 2 JJ|^| V _ 3 < C(V, d, 1/)) E\g \ v ]J \gj\v-i- 

3=3 3=1 

The estimates given above show that (|3.15|) holds for v = V as soon as it holds for v = V — 2 and 
v = V — 1. We have already proved that (|3.15|) holds for v = 0,1. Thus, (|3.15|) holds for every v. 
The lemma is proved. 

Proof of Lemma We have already proved (|3.10|) to hold for every g G Q H . Now, let g G 
W^+^iH). Consider {g n } C G H such that g n -» g in W^ l ~\H) (recall that Q H is dense in any 
Wf(H) by definition). By (IBTTUD . for any k = 0,...,l, 

\\D k 5(g n ) - D k 5(g N )\\ L2m ^ P:H m) < C(m, I, d, ip)\\g n - 9n hm,2m+i-i 0, n,N ^ +00. 
Thus there exist F k £ L 2m (VL, P, H m ), k = Q,...,l such that 

\\D k S(g n ) - F k \\ L2m{Qi ^ Hm) -> 0, n-^+oo, fc = 0, . . . , Z. 
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Since operator 5 is closed, Fq = 5(g). Using that operator D is closed, one can verify inductively 
that F k = DF h _i, k = l,...,l. This means that 6(g) G W l 2m with D k 5(g) = F k , k = 0, . . . , I. At last, 
using (|3.1(jp we get 

I i 
ll<%)fc = H^fc ^ hmsupi^ \\D k 5(g n )\\ 2 ™ m = 

k=0 n k=0 

= limsup||% n )||l^ < C(m,l,d,i/})Um sup ||#n||im,3m-l = C (™> d i ^hWlm^m+i-l- 
n n 

The lemma is proved. 

3.3. Malliavin's representation for the densities of the truncated distributions of smooth 
functionals. The typical result in the Malliavin calculus on the Wiener space is that, when the 
components fi, ■ ■ ■ , fd of a random vector / = (/i, . . . , /^) are smooth enough and the Malliavin 
matrix a* = {(Dfi,Dfj)n}fj =l is non-degenerate in a sense that 

(3.21) [det^p 1 G p| L p (Sl,$,P), 

the distribution of / has a smooth density (see, for instance, [12] §3.2). Such kind of a result is useless 
in the framework, introduced in subsection 3.1, since there does not exist any functional / satisfying 
(|3.2ip : if s\ = ■ ■ ■ = e n = then Df = for every / G C. In order to overcome this difficulty we 
use the following truncation procedure: we consider, instead of P, a new (non-probability) measure 
Pz(-) = P(- H H) with some set 5 G a(e, C). If this set is chosen in such a way that (|3.2ip holds true 
with P replaced by P~ then the Malliavin's calculus can be applied in order to investigate the law 
of / w.r.t. -P-. In this subsection, we give the Malliavin's representation for the density of this law. 
All principal steps in our consideration are analogous to those in the standard Malliavin calculus on 
the Wiener space (see, for instance, [12] , Chapter 3). Therefore, we sketch the proofs only. 
Let fx, . . . , fa G C be fixed, consider the Malliavin matrix a? = (a{j) d j =1 , 

a{j = (Df t ,Dfj) H = X>fafc)[Qfer/i(fae,C))] ■ 0kr/i(fa,e,C))]- 

k.r 

Consider a set H G a(e, () such that E C {det a* > 0} and 

(3.22) El s [deta f ]~ p < oo, p > 1. 
Put 

' [a- f (lu)}- 1 , wGS 

0, LU(?Z' 



gf> s (u) 



Proposition 3.4. gf' S G W™(R dxd ) and 

(3.23) (Dg f > s ,h) H = -g f > B ■ (Da f ,h) H • g f ' B , heH. 

Sketch of the proof. It is enough to prove that gf' 3 G f] P >iW p 1 (R dxd ) and (l3~23l) holds true. 
Suppose that o~* > d K d with some c > 0. Then one can easily see that g^^ G C and (|3.23p follows 
from the well known formula for the derivative of the inverse matrix, 

j t [A(t)]- 1 = -[A(t)]- 1 • [j t A(t)\ ■ [A(t)}-\ 

In the general case, consider the matrix- valued functions a^ c = a* + c/ R d and g^~-> c = 1= • [cr^' c ] _1 , 
c > 0. Condition (l3~22l) provides that g f ^' c g f ' E ,c -> 0+ in any L p . It is already proved that 
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(|3.23p holds true for the functionals indexed by c. Thus, passing to the limit as c — > 0+, we obtain 
the required statement. 

Denote $p = Ylt=i @kT ' Dfki i = 1, ■ ■ ■ ,d. Also denote, by Eg, the expectation w.r.t. P=. 



Proposition 3.5. For every i = 1, . . . , d, u£ W££ and every F G 
(3.24) E*[diF]{f u ...,f d )-v = E s F(h, . . . , f d ) ■ S (v • tfp 



Sketch of the proof. It follows from Propositions l3.3tl3.4l and Lemma lBTTl that v-fif G Dom{5). Since 
S G a(e, £), the function 1= belongs to C and has its stochastic derivative equal to 0. Proposition 13.21 
provides that <5(ls • g) = l^(g),g G Dom(5). Therefore 

d 

E[diF]{f u ...,f d )-l B -v = J2 EfoFKh, . . . , f d ) ■ I s • v ■ [at ■ = 

3=1 



d d 

J.S / 

j=l k=l 



EE^M(/i,...,/d)-Ia-«-fifi a -oj 



d - / d 

EJ2 Is • v ■ g(fi Y^[djF]Ui, ■ ■ ■ , /d) ■ D fj, D fk = E (e>[F(/i, . . . , f d )],l s ■ v ■ 
fc=i \j=i / H 



ii 



= EF(/ 1 ,...,/ (J )- ( 5(l 3 -i,-tff) =£F(/ 1 ,...,/ d )-l H ^(^^), 

that provides (pHl")) . 
Put 

v('~ = 1, = S (y{' s • tfp) , I = 1, . . . , d, 

T f '* = vg v Tp = 6 (V> 3 • tfp) , i = l,...,A 
Write P~ for the distribution of / w.r.t. P«. For ai, . . . , a d G {0, 1}, denote 

lai...a d (x) = l(-l) Q ixi>0,...,(-l) a d Xd >0) x G ^ 

Proposition 3.6. The distribution P~ /ias a density pt, bounded together with all its derivatives 

f 

dip s ,i = 1, . . . ,d. For any ai, . . . , a d G {0, 1}, 

(3.25) p /( y ) = (_i)-i+-+« djE ; lQi ^ ad ( / _y).T/' s , yGM d , 

(3.26) d iP f s (y) = (-ir + - +ad+a * +1 El ai ... ad (f - y) • Tp, y G R d . 

Sketch of the proof. Applying iteratively (|3.24p one can deduce that, for every F G C£°(M d ), 

(3.27) Esfo . . . d d F](f 1} . . . , f d ) = EF(fi, . . . , f d ) • T^> s , 

(3.28) E s [dt . . . d d diF]{f u ...,f d ) = EF(fi, . . . , f d ) • Tp, i = 1, . . . ,d. 
Now, the informal way to get representation (|3.25 j) is to apply (|3.27p to F = l ai ... ad - 

p f 3 (y) = . . . d d E B l ai ... ad (f -y) = (-ir + -^E E [d 1 . . . d d l ai ... ad ](f - y) = 

(3.29) = {-ir + - +ad El ai ... ai {f - y) ■ 
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In order to justify (|3.29p one should consider smooth approximations F n for the function F = l Ql ... Qd 
and use Fubini theorem (we omit detailed exposition here, referring the reader, for instance, to |12j . 
§§3.1, 3.2). Similarly, (|3.26p is provided by (|3.28|) and the formula 

d iP l(y) = (-ir + -+ Qd+1 £ H [di . . . d d dd ai ... ad ]{f - V). 

3.4. Estimates for the densities of the truncated distributions of smooth functionals. 

Proposition 13.61 immediately provides the following family of estimates for the density p~ of the 
truncated distribution of /. 

Corollary 3.1. For any y £ M. d , 

pl{y) < ||T^ S || L2 • min pj ((-ir/i > • • • , (-i) ad fd > (-i) ad y d ) < \\T f ' s \\ L2 , 



tL{v) < \\T f rh 2 • min P| > (-l) ai yi, (-l) ad f d > (-l) ad y d ) < ||Tp| 



i = 1, . . . ,d. In particular, pL satisfies Lipschitz condition with the constant L = X]j=i ll^i \\l 2 - 

In this subsection we give explicit estimates for ||T^ ,,= '||i 2 , ||T^'"[|x 2 , i = l,...,d. Our estimates 
somewhat differ from the standard Malliavin-type ones. In our considerations, we operate with 
the matrix [cr-^] -1 straightforwardly and do not use (unlike in the standard Maliavin's approach) 
representation of this matrix via the Cramer's formula [cr-^] -1 = [det cr^] _1 • S-^ (S-^ denotes the 
cofactor matrix for a* ). This is caused by our goal to prove, together with existence of the density, 
an explicit estimates for it like the estimate (ii) of Theorem 11.11 

Let us give an iterative description of the family {vf' a } involved into construction of T^ ,a ,T^' a . 
We introduce two families of operators acting on W^: 

Ii : ip i-> 5(ipDfi), J ijk : ip ip(Da f jk ,Dfi) H , i,j,k = l,...,d. 

We call any operator 1%, . . . , I d an operator of the type I, and any operator from the set {Jijk, i,j, k = 
1, . . . , d} an operator of the type J. We denote by JC(m, M) the class of all functions that can be 
obtained from (p = 1 by applying, in arbitrary order, of m operators of the type / and M operators 
of the type J. 

Proposition 3.7. For any I = 1, . . . , d + 1, there exist constant C(d, !) 6ff such that is a sum 
of at most C(d,l) summands of the type 



(3.30) <p Jl g 



r 

u 

k=l 



where ik,jk = 1, ■ ■ ■ ,d are arbitrary and ip belongs to some class 3C(m, M) with in + M = I — 1 and 
r = M + 1-1. 

Proof. We use induction by I. For / = 1, the statement is trivial since v(' a = 1 6 %(0, 0). Suppose 
the statement of the Lemma to hold true for some / < d. Let us prove this statement for I + 1. Due 
to the inductive supposition, is a sum of at most C(d, I) summands of the type 



k=l 



20 ALEXEY M. KULIK 

ip G X(m, M),m + M = l,M + I = r. We have 

s {v ■ fl = E <^ • [fl eSj • # • ^ 

fe=l <J=1 fc=l 

Thus, Uj is a sum of at most d ■ C(d, I) summands of the type 

r+l 



k=l 



Lp E X(m,M),m + M = l,M + 1 = r. Due to Propositions E21 and E 



r+l 



r+l 



r+l 



fc=i 



/c=i 



r+l 



l+l 



(3.31) 



<?=! 



k=l 



n e 



/.B 



k=l 



gf<* -(Dat,Dfi)- g 



H 



f 



k<r+l,k=£q 

Since (p E X(m, M), 5{<p> ■ Dfi) G 3C(m + 1, M). Thus, the first term in the right hand side of (I3.31D 
has the form (|3.30p . Every summand in the sum in the right hand side of (|3.3ip is a sum of d 2 terms 
of the type 

>+2 



ip-(.D<4j,Df,) H 



fc=i 
/ 



ikjk 



Every such a term has the form (I3.30|) . since <p ■ (Dai-, Dfi)n G X(m,M + 1) for ip E X(m,M). 

ij 

Therefore, the statement of the Lemma holds true for l + l, also, with C(d, l+l) = C(d, l)[l + d 2 (l + l)]. 
The proposition is proved. 

Recall that Y^'~ = vf ,ad+1 , and thus Proposition 13.71 provides that T^' a is a sum of not more than 
C(d) summands of the type (|3.30|) with <p> E %{d — M, M) and r = M + d (M may vary from to d). 
For every such a summand, 

M+d 

(3.32) 



II &U 2 < IMu 



^WP \\M 



k=l 



where ||^4||a4 = maxjj^,..^ \Aij\, A E M. dxd . Thus, in order to estimate ||T^' ,= '||i 2 , it is sufficient to 
estimate rnax^xid-MM) IMU 4 - Denote on = Dfi,[5 ijk = (Da f jk , Dfi) H . 



d-M 



Proposition 3.8. For every ip E %{d — M, M), 
(3.33) \\<p\\u < • ( max||a i || 2 (d+i)( ( i+2),(d+i)2-i 

Proof. By Lemma [3TT1 and Proposition 13.31 f° r an Y * = 1? • • • j d, (p G W 7 "^, m > 



max ||/%fc|| 2 (<f+l)(<i+2),(d+l) 2 -l 

ijk 



\\${V a i)hm+2,m?-l < d , ip) \\ipaa || 2 m+2,(m+l)2- 



< 



(3.34) < C(m, d, V , )ll9 9 ll2m+4,(m+l) 2 -lll a «ll2(m+l)(m+2),(m+l) 2 -l- 

By Proposition 13.31 for any i,j,k=l,...,d,(p£ W££, m > 



\\ i fPijk\\2m+2,m 2 -l — C(m) |Ml2ra+4,m 2 -lA.jfc||2(m+l)(m+2),m 2 -l — 
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I 3 - 35 ) < C(m, (i, V')||9 ::> ll2m+4,(m+l) 2 -ll|Aifc||2(m+l)(m+2),(m+l) 2 -l- 

Recall that || • = \\ ■ 1 1 an d £ %{d — M,M) is obtained from 1 by applying (in some order) of 
d — M operators of the type I and M operators of the type J. Thus, in order to obtain (|3.33j) . one 
should consequently put m = d, d — 1 , . . . , 1 and apply either inequality (j3.34j) or inequality (13,351) 
depending on what type of the operator (/ or J) was applied at this position in the construction of 
the function ip. The proposition is proved. 
J 



Recall that aj k = {Dfj,Dfk)H- By Proposition 13.3 



\Pijk\m < C(m)\Dfi\ m \Da f jk \ m < C(m)\Dfi\ m \Dfj\ m+ x\Df k \ m+ i. 



Thus, Proposition 13.81 provides the following estimate for ||T^ ,= ||i 2 . Denote 
N d (f) 



max max 
i=l,...,dm=l,...,(d+l) 2 L 



E\\D m h 



,2(d+l)(d+2) 



2(d+l)(d+2) 



Corollary 3.2. There exists a constant L d , dependent on d and if) only, such that 

d 



(3.36) 



||T^|| L2 < l d £ [N d (fj 

M=0 



d+2M 



^we \\m 



For ||T^'"||i 2 , the following estimate holds true (we omit the proof since it is totally analogous to 
the proof of (|3.36p given above). 



Proposition 3.9. For any i = 1, . . . , d, 

d+l 

\ L2 <i d+ i £ [;v*fi(/) 

M=0 



(3.37) 



IT 



d+2M+l 



E\\g 



/,S||4(M+d+l) 



At the end of this section, we formulate a general local limit theorem. This theorem is a straight- 
forward corollary of the representation given by Proposition 13.61 and the estimates (|3.36p and (13.370 . 
Consider the sequence of probability spaces {(f2 n , SF 1 , P n ), n > 1} of the type (|3.3|) with the given 
measures ttjj, x, v. Let the function if) and the functions n be fixed. Denote H n = W dxn and consider 
the derivative, gradient and Sobolev spaces constructed in subsection 3.1. For a sequence of random 
vectors f n :Q n ^- M. d and a sequence of sets {H n G a(s, Q} denote by P n ^ n the distribution of f n 
w.r.t. P n and by P« the distribution of f n w.r.t. P£ = P n {- n H n ). Denote 

d 



K d (/,S) 



d+2M 



E\\p f;L 



MM+d) 



M=0 



Theorem 3.1. Suppose that {f n } and {S n } satisfy condition 
(CI) f n € W^ {d+3) (M d ) and sup n K d+1 (r,E n ) < +00. 



Then P, 



possess 



a densities p~ . Moreover, 



(a) pi n (y)<K d (r,E n ).Pl(\\r\\>\\y\\); 

(b) p4 n satisfy Lipschitz condition with the common constant equal to d • sup n K d+ \(f n , S n ). 
//, additionally, 

(C2) f n converge in distribution to some random vector f; 
(C3) P n (~„)^l,rw+cx>; 
then the distribution of the vector f possess a density pf and 
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(c) sup^jjd \pt (y ) - P f (y ) | -» 0, n -> +00 . 



Proof. Although the statements of Corollaries 13.11 and 13.21 are formulated for / £ C Rd , they can 
be extended to / £ W^ftf) ) ^ f ,,^ ) (M d ) by a standard approximation procedure. For any y £ R rf , there 



2(d+2)(d+3) v 

exist a choice of the signs a±, . . . ,ad such that 

ph„ ((-ir/r > (-ir^i, . . . , (-lr-/? > (-ir d y d ) < Ps n (\\n\ > my 

Therefore, statements (a) and (b) follow immediately from Corollaries 13.11 and 13.21 Statement (b) 
provides that the sequence {pt n } has a compact closure in the space C(M rf ) with the topology of 
uniform convergence on a compacts. This together with the conditions (C2),(C3) provides (c). The 
theorem is proved. 

4. Proofs of Theorems II. II -[2~2l 



4.1. Proof of Theorem 11.11 We reduce the proof of Theorem ll.ll to the verification of the conditions 
of Theorem 13 . 1 1 and explicit estimation of the expression in the right hand side of (a). We use, without 
additional discussion, notation introduced in Section 3. 

Denote f™ t = X n (t) — x, where the processes X n are defined by (|0.2p . (IO,3p with the initial value 
X n (0) = x £ When it does not cause misunderstanding, we omit the indices x,t and write f n 
for fl t . 

We conduct the proof in several steps. First, we give explicit expressions for the derivatives of 
the functionals f n . Next, we estimate the moments of these derivatives (this allows us to estimate 
Nd+i(f n ))- Then, on the properly chosen S n , we estimate the inverse matrix for the Malliavin 
matrix o~f™ (this allows us to estimate K^_|_i(/ n , H n )). At last, we estimate the tail probabilities 
-ffe n (||/ n || > \\y\\) m order to provide the estimates given in the statement (ii) of Theorem ll.il 

Everywhere below we suppose conditions (Bl), (B2 K ), (B3) of Theorem ll.ll to hold true. We prove 
(i) - (hi) in details and give a brief sketch of changes that should be made in order to prove (h'),(iii / ). 
We put 

K — 1 

(4.1) n (C) = I max fc < n ||Cfe||<n«, ? = 2k + 2 ' 

In order to make notation more convenient, we rewrite (10. 2D to the form 



X n 

n I \ n I \ \ n J J n z — ' \ \ n I I \ln 

' x ' x x ' ' r=l v v / / y 

here . . . , £kd are the components of the vector £/% and b±,...,bd are the columns of the matrix b. 



Lemma 4.1. For every t £ [0, 1], X n {t) £ f| p >i Wp ( M )• Derivatives Y n (t) = DX n (t),t€ [0, 1] 
satisfy relations 

K„ (0 ) _ 0, r n (5) _ y. + v. (x n (^i)) y. . U 



(«) + ± [v, (¥)) * ( V) ■% + (,„ (izi) 



® ekr 
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k = 1, . . . , n. 



Sketch of the proof. The proof is quite standard, and thus we just outline its main steps. Using 
induction by k, one can easily verify that, for every j,k,r, there exists dj r X n (^) = (Y n (^) ,^jr)H n 
with Y n defined by A4.3H . One can see that Y n = as soon as maxfc< n ||Cfc|| > n q and, there- 
fore, ess sup || Y n (~) || < +oo for every k < n. Iterating these considerations, one can verify 
that ess sup ||V^X n (^) || < +oo for every k < n,m < (d + 2) 2 , that means that X n (^) G 



n P >i^ 



(d+2) 2 



with DX n (&\ = ^n(^)i that gives the statement of the Lemma for t 



For arbitrary t £ [0, 1], this statement holds by linearity. 
Denote = E\\& \\ K . 

Lemma 4.2. For every p > l,m£N, there exist constant C(a,b,d,U, fi K (^),m,p) such that 
(4.5) E\\D m X n (t)f m ^ t <C( a ,b,d,U,^(0,m,p)-ti, !€ 



Proof. Consider first the case m = 1. It is enough to prove (14. 5h for t = = l,...,n and 
p = 2q, q £ N. We have 



Y, 



k-1 



+ 1 -(H— 

n \ V n 



k-1 



, Va X, 



- 1 



fc - 1 



n 



1 



+ ^(Va(x ri 

r=l ^ ^ 



n,r2=l 



n 

n 

k-1 

n 



Va X 



Yn 



Y n 



Y, 



k-1 



71 



k-1 



n 



,V6 r X 



, V6 r2 X n 



A; — 1 



n 



n 



k-1 



Y, 



Yn 



+ 

k-1 

n 

k-1 

n 



Y n 

/n 



fc- 1 



n 



+ 



+ 



+ 



n 



(4.6) 



+ 



^(C)i £fc =iV(%) 



6, X 



k-1 



n 



here we have used the fact that Y n , Va (X n (*=±)) F n , V6 r (X n (^)) Y r 

1, . . . , d belong to the subspace generated by the vectors of the type v®ej r ,v G M. d , j < k,r = 1, . . . , d 
and b r (X n (^p)) &> efc r , r = 1, . . . , d are orthogonal to this subspace. 

Denote \\Y n (^) H^m^ = T it- Recall that Y n (^) = 0, /c = 1, . . . , n as soon as there exist j = 1, . . . , d 
such that ||Cj || > re 11 . Since coefficients a, 6 are bounded together with their derivatives, we can rewrite 
gSD as 



1 V- &r V- &n&r 2 

T fe _i + 6 fc _i • - + 2^ A *-V " 77^ ' I IICfcll<n' ; + A fc-i,n,r 2 1 

r=l v ri,r2=l 



IICfell<n 5 



xl 



/c = 1, . . . , n 
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with S'j fc _ 1 = a(rn,ei,Ci, ■ ■ ■ ,Vk-i,£k-i,Ck-i) - measurable 6 fc _i, Afc-i,n,r 2 such that 
(4.7) |6 fc _i| <C(a,6,d)[l + T & _ 1 ], |A fc _ liPlira | < C(a,6,d)T fe _ 1 . 

Since Y*. > 0, we have 



(4.8) 



i=0 



n 



' A IICkll<n< 



We have ^ r = r/fc r + ^ r and the set C7 of the possible values of r/fc = (?7fci, . . . , rjkd) is bounded. In 
addition, |Cfcr|I|[c fc |l<n« — n ' '• Therefore, 



(4.9) 



fcr 



" I IIC fe ||<^ < 



(C(U) ■ n 

I 

J12 



C l (U)n l ^)<^l, l> K + l 



Since £ ; ||Cfc||' t < +°o, 
(4.10) E 
Similarly, 
(4.11) 
At last, 



5 kr 



■ I IIC fe ||<^ ^ 



^ £(||£ fe fVl) 



< 



n 



l = 2,...,K. 



E 



n 



l IICfell<^ 



< -^II6|| 2 < 

n 



£(||&|| K vi) 



n 



E 



n 



l IICfell<™^ 



E 



£,kr 



71 



l IICfell>^ 



(4.12) 



<C(U,fj, K (0)-n—* 



< n 2 



n 



E\^kr\ K \ K [P{Kkr\\ > n ? ) 



< 



< 



n 



(recall that E^r = 0). The triple Cfc) is independent of S'k-i- Thus, taking in (|4.ip conditional 

expectation w.r.t. 3~fc_i and taking into account inequalities (|4.7f) . we obtain an estimate 



-ET? < E ( T fc _! + 9*_i • - ) +C(a,b,d,U,ii K (£),q) 



n 



n 



< 1 + 



Ci(a,b,d,U,n K (£),q) 



■n 



(4.13) 



9-1 / 1\ 

+ C 2 (a,6,d,?7,^(e),Q)E( 1 + -) 

z=o ^ n ' 



IV 1 



ni- 



Let us show that (|4.13p provide the family of estimates 

'k^ q 



(4.14) 



Erl<C(a,b,d,U,fM K (0,q) 



n 



k = 1, ... ,n, g£N 



(note that (|4.14p is exactly (14. 5|) with m = 1 and p = 2q). We use induction by q. For g = 1, (I4.13|) 
implies that 

Ci(a,6, ci,C/,// K (£),!)' 



^ , C 2 (a,b,d,U,n K (£),l) C 2 (a,b,d,U,n K (£),l) 
& 1 fc S 1 



n 



1 + 



+ ••• + 
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C 2 (a,b,d,U,n K (£,),l) f Ci(a,6, d, U,/j, k (£), 1) 



k-l 



n \ n J n 

Similarly, if (|4. 14|) holds true for all q < Q — 1, then (|4.13p implies that 



< -■C 2 (a,b,d,U,v K (0A)-e Cl{a ' b ' d ' U '^ m) - 



£T« < C 2 (a,b,d,U,» K (Z),Q) ■ e W,dMO,Q) . £ . C (a,M,£W£),0 • (~ ) < 



z=o 



< c(a,5,d,[/,^(e),Q) 

that proves ()4.14p for q = Q. This proves the statement of the lemma for m = 1. For arbitrary ra, the 
proof is analogous: one should write difference relations for the higher derivatives of X n , analogous 
to (|4.3p . and then again use the moment estimates of the same type with the given above. This step 
does not differ principally from the one for SDE's driven by a Wiener process (see, for instance |13j . 
Chapter V §8), and thus we omit its detailed exposition here. The lemma is proved. 



Corollary 4.1. 



Nj(X n {t)) < C(a,b,d,U,n K (Z),j) -Vi, t€ 



n 



,jeN. 



Let us proceed with the investigation of the Malliavin's matrix cr-^" for f n = X n (t). Denote by 
£™j 5 < i < j < n the difference analogue of the stochastic exponent for (14. 3p . i.e., the family of 
M. dxd - valued variables satisfying the relations 
(4.15) 



j-1 



n 



1 d 



r=l 



V/), ( X„ ( - — -Y) 



n 



j = i, . . . ,n. Then one can easily obtain the representation for Y n (-) 

k d 



(4.16) 



Yn 



j = l T = l V 



j-1 



&kr- 



Denote f n > k = X (J) and a nyk = a fn ' k . By (fTOjh we have 



k d 



j=l r =l 
k 



E% k br ( X 



n 



^[^(^•)l £ ,=i] 



££ fc 6 ( X n 



J-1 



n 



.(0 



Y,[^(Vk)le k =l]-El k -[bb*](x n 



3=1 



j-1 



^j,k^r ( X n 



J-1 



7? 



J-1 



71 



77 



Together with the family {£",}, we consider the family defined by 

(4.17) 



r=l 



V6 r X 



J-1 



77 
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j = i, . . . , n, where £™ = ^l||^||< n «,i = 1, . . . ,n. By the construction, £"^- = on the set {6 n (() 
1}, therefore 

= ^ ^[v 2 (^)i £j =i] • q fc • [»*] (x n ( 3 —)) • [£> fc ]*. 



n 



3=1 



11 



We have 



i+l 
II 



'•J 



+ v "( x -( i i i ))4 + g vi "( x " 



I - 1 



SZr_ 



Since Va, Vo are bounded and |£™ r |I||f -||<n« — max(max xg [/ ||x||, n q ), there exists no = no(a, b, d, U, <; 
such that 

d 



Va X 



,_1 )H+E^(*. 

/ / r= i \ 



71 



I||£||<n* < ^, n > n . 



Then £f • is invertible and 



/rod + Va X 



fer 1 = n 

l=i+l 

Thus, on the set {6 n (() = l}\{ei 



?i 



n / / v/n 



£n = 0}, the matrix a n ^ is invertible and 



(4.18) 



n,k I 



inf K^^)]- 1 < [max life]" 1 !!] 2 • J i^^ 2 fe)I £ . =1 1 



Lemma 4.3. For every p > 1, 



£[max || [e^]" 1 !!^ <C(a,6,d,C/,^(0>P), n € N. 



i<j<n 



Proof. Since [£"•] 1 = £q^ • [£o j] > it is enough to prove that 
(4.19) S[max||£^||F < C(a,b,d,U, » K (0,p), 



i<n 



(4.20) 

Let us prove inequality 
(4.21) 



.E max | det £q ■ | p < C(a,b,d,U, /j, K (£),p). 

j<n 



E[max||£^|| 2 f < C(a,b,d,U,fji K (0,p) 

i<n 



with \\A\\ 2 = yj A tr'> this wil1 provide inequality (|4TT§j) . We deduce from KTt} that Zf = ||£^ 
satisfy relations analogous to (|4.6p . i.e. 



(4.22) 



1 £ n £ n £ n 



ri ,ra 



with an {3~j} - adapted sequences V^' n , V^ 2 ' n , V^ 3,ra such that 

(4.23) \Vl> n \<C(a,b,d,U)(l + Z?), \V?; n \ < C(a, 6, d, 17)2? , |^ r J < C(a, 6, d, U)Z? . 
Then then the moment estimates analogous to those made in the proof of Lemma 14.21 provide that 



(4.24) 



max(EZ?)* <C(a,b,d,U,ii K {Z),p). 
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3 A An A/Tn _ Y^i 



Denote A* = Z% + YX=x AA i > M j™ = £*=l AM ™ with 

era TPC n 



and 



1 ~FF n Ff n f n 

r » 

Then Z< = Aj + M<. By gT) and (Il23| . 



| AA n|<C(g,6,^^(0) -(1 + zo> 



and therefore (|4.24|) provides that 



Emax.\A5\2 < C(a,b,U, fj, K (£),p). 

Similarly, Burkholder inequality together with (|4.23|) and (|4.9p - ()4. 1[) provides that 

Emax\M?\% < C(a,b,U, fi K ^),p), 

j<n 

that proves (|4.21|) . and therefore (|4.19j) . 

On the set {||A||2 < \} C R dxd , the function $:4h (det[/ R d + A])' 1 can be represented in the 
form 

$(A) = 1 + Q(j1)+0(A), 
where P is a polynomial of A with degA < nd and |$(A)| < C||A||2 +1 - We have 

C(a,b,d,U) 



Va X. 



l-l 



n 



)H+x>H^))-§ 



< 



Therefore, 



(detfi^r^CdetE^.!)- 1 



i + Ql 



era era 
sil sjd 



n' ' \/n 



where |#"| < ; #«■ j s g^. _ measurable, Q?_i is a polynomial with degQj-i < and its 

coefficients are 3^-1 - measurable and bounded by some constant depending on the coefficients a, b. 
Repeating the arguments used in the proof of (|4.19p we obtain (I4,20|) , The lemma is proved. 

Lemma 4.4. For every p 6 N, c > 0, 
k 



E I 



< C(c,V',p)fc~ p , k> 



2p + l 



Remark. For arbitrary ip £ C°°(IR rf ) with i/j = on <9J7, the given above statement may fail. It is 
crucial for ifi to have non-zero normal derivative at (some part of) the boundary in order to provide 

Proof. Since r\ and e are independent, 



}ck\ 



j'=i 
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where ]x]f= min{n G Z|n > x}. By the construction of the function ip, 

(4.25) P{ip 2 {r}j) < z) ~ C(V>) • v/I, z ->■ + . 
Therefore 

P(4>\m) + ■■■ + if?(w) < z) ~ • z o+ 

and 

(4.26) ^ 2 (r?i) + --- + ^ 2 (r/0) _P < +oo 



as soon as I > 2p. We put q = 2p + 1 , N 



]ck\ 



and divide the set {1, . . . , ]ck[} on the blocks 



{1, ...,q}, {q + 1, ... , 2q}, . . . , {(N - l)q + 1, . . . , JVg}, {TVg + 1, . . . ,}ck[} 
(the last block may be empty). We denote 

iq 

Y, V> 2 (^), i = l,...,N. 



j=(t-l)g+l 



We have 



f ]<*[ ] p / N \~p / N \ 



The function x i— > x p is convex on M + , and therefore 

-p 



£??]~ p < +oo 



(the last inequality follows from (|4.26p ). If k > ^j" 1 , then > 1 and therefore N 



}ck[ 



4p+2 



fc. Thus 



MI 
<? 



^ 2q ^ 



}ck[ 



4p + 2 



The lemma is proved. 

Inequality (|4.18p . Lemma 14.31 and Lemma 14.41 provide the following estimate. For a given c > 
and t G [0, 1], we put S n = {# n (0 = 1}D {eJS e i > c N}- 

Corollary 4.2. Forp G N and [in] > 2£ ±i, 



Bp' 



<C(a,b,d,U^,^(0,p)-t- p . 



At last, let us give an estimates for the tail probabilities for f n . The following lemma is completely 
analogous to Lemma 14.21 the proof is omitted. 



Lemma 4.5. For every p > 1, i/iere exist constant C(a,b,d,U, p, K (£),p) such that 
E\\X n (t) - X(0)||; a I,,. (c)=1) < C(a,b,d,U,p,Ai),P) • * ! , * 6 



n 



Corollary 4.3. For every p > 1, i/tere exists constant C p , dependent on a,b,d,U, n K (£),P> such that 



p(\\x n (t) - x n (o)|| > i/,e n (0 = i) < c p i + 



2\ -P 



n G N,£ G 



n 
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Lemma 4.6. There exist constants Ci,C2,C3 ; dependent on a,b,U,d, // K (£), such that, for every 
A G R d with ||A|| < Cin 7 ^ 1 , 

(4.27) Ee^ x ^- X ^h {en{0=1} < C 2 e c ^\ t G 



1 

- 1 

n 



,n G N. 



Proo/. For a given A, denote Z n (t) = e^ X ' Xn ^ We have Z n (0) = 1. On the other hand, 



■shr 



< 



m&x yeU \\y\ 



+ n ■ n 2 



max yeC / ||y| 



on the set {#n(C) = !}• Thus there exists a constant C4 such that, for ||A|| < Cin re+1 . 

d 



- ( A. </ j A',., 



fc-1 



n 



k-1 



n 



<C 4 



on the set {9 n (0 = !}• Using the elementary inequality e x < 1 + x + Cx 2 , |x| < C4, we obtain that, 
on the same set, 

d 



Z„ [ - 

n 



Zn 



< z„ 



k-1 



n 



cxp 



— I A. (/ ( A,, 

n 



k-1 



n 



k-1 



1 + &k-i 



|A|| + IIAI 



n 



(4.28) 



~l~ ^ ^k—l,ri,T2 
ri,ra=l 



d 
r=l 

£fcri ^fcr2 1 1 ^ 1 1 



k-1 



n 



< 



MI|A|| + ||A|| 2 ) 
fe-l,r " /= 1" 

n 



?i 



with an S^-i - measurable coefficients Ofc_i, Ajt_i n r2 , bounded by some constant C. An arguments, 
analogous to those used in the proof of Lemma 14.21 provide that (|4.28p implies the estimate 

|2\ k 



(4.29) 



EZ„ [ \ ) < [ J + C 



A + 



< exp 



2C .-.(1 + A 2 ) 

n 



,k,n G N. 



This is exactly (I4.27|) for t = |. For i G (f ,^), Zn(*) is a linear combination of Z n (£) and 
Z n (^jp). Therefore, (|Q7j) follows from ([09]) via Jensen's inequality and relation < 2% < 2t 
(recall that t > — and thus k > 1). The lemma is proved. 

Corollary 4.4. There exist constants C^,Qq,Qt , dependent on a,b,U,d, swc/i i/iai 



and 



P(\\X n (t) - X n (0)\\ > y,9 n (C) = 1) < C 6 e- C7 ^, y G (0,C 5 in^),n G N,t G 



P(||X„(t) - X n (0)|| > j,,0 n (O = 1) < C 6 e- C7n ^ y , y > C 5 in^,n G N,t G 



1 



n 



n 



-,1 



Proof. It is enough to verify that, for any coordinate (X n )j of the process X n , j = 1, . . . ,d, there 
exist constants C5, C6, C7, Cs such that 



(4.30) P(±((X n ) j (t) - (X n )j(0)) > y,9 n (C) = 1) < C 6 e- C7 ^, y G (0,C 5 tn^),n G N,t G 



1 

-,1 

n 
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and 

(4.31) P(±((X n ) J (t) - {X n )j{0)) > y,e n (C) = 1) < Qe-^^y, y > C 5 to^,nG N,t € 



1 

- 1 

n 



Inequality (^"30]) with C 5 = 2CiC 3 and C 7 = [2C 3 ] _1 follows from (jOT)) with A = (±2^] -e^, where 
ej is the j-th coordinate vector in R d . Inequality (|4.3ip with the same C5 and Q = % follows from 
(gSZD with A = ('rbCm^+r) • ey. 

Proof of Theorem ] l.li We take p = 8(d + 1) and fix some c S (0, a) (a is given in condition (B3)). 
We write n* = no (a, b, U, (see the notation before Lemma |4.3|) and put 

R(C) = 1} f1 {Eg £j > c[tn]} , n > n., [in] > 

otherwise 



QlMv) = P(X n (x,t) G dy,Ei) = P &n (fl t edy-x), 
Rlt(dy) = P(X n (x,t) e dy,n\H^). 

Corollaries 14.11 and 14,21 provide condition (CI) of the Theorem 13.11 By the statement (a) of this 
theorem, 

(4.32) Ql t (dy) = qlMdy with g« t < K d (f n ,E n ) • P^{\\f n \\ > \\y\\,9 n (0 = 1). 

Moreover, Corollaries 14.11 and 14.21 provide an explicit estimate for K^(/ n ,H n ). Namely, for some 
constant C dependent on a, b, c, d, « K (£), U, ijj, 



(4.33) Kd(r,Sn) < C J] [^] d+2M • [£~ 4 ^+ d )]i = (d + l)Ct"i 

M=0 

Thus the statement (a) of Theorem 13.11 and Corollaries 14.414.31 provide statement (ii) of Theorem 11.11 
By Chebyshev inequality, 

P(6 n (() = 0) < n ■ n-^ = n~ e{K) . 
Take A = In ( °^Za) ) - > 0- By Chebyshev inequality, we get, after some simple calculations, 

(4.34) P ]J2 £ i < ck J < -^xck = [*(«' c )] > k G N ' 

with *(a, c) = (^) C - One can verify that ^(a, c) < for < c < a < 1. Thus, we can 

conclude that, for p = — \ In ^(a, c) > 0, 

P(0\5*) < n" e(K) + e- pn * when n > n*, [tn] > ?2±1 

c 

(we have used here that [nt] > ^ for t > -). Thus, for all £ > 0, 

(4.35) ^(0\3^) < £>[n~ e(K) + e"""*] 

with the constant D dependent on n* , p, c, p. This provides statement (iii) of Theorem 11.11 

We have shown that if x £ M. d , t > are fixed then the functions f n = f£ t and the sets E n = 
satisfy all the conditions of Theorem 13. 11 This means that (Jxtiv) ~^ Px,t(y) uniformly w.r.t. y 6 M. d . 
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In order to show that this convergence holds uniformly w.r.t. t G [6, l],x, y G M. d , we need to show 
that, for every sequences {t n } C [6, 1], {x n }, {y n } C R d , 

( 4 -36) q n XnU (y n ) ~ Px n ,t n (y n ) -► 0, n -> oo. 

We can suppose that {t n } converges to some i G [5, 1]. The functions a, 6 are bounded together with 
their derivatives up to the second order, and therefore the sequences of the functions 

a n (-) = a(- + x n ), b n (-) = b(- + x n ) 

are pre-compact in C 1 (M, d , and C 1 (R d , M rfxrf ), correspondingly. We can suppose that 

a n -» a in C^R^R"*), b n -» b in C 1 ^, M dxd ). 

Consider the processes Z n defined by the relations of the type (|0.2p . (|0.3p with Z n (0) = and the 
coefficients a, b replaced by a n , b n . Also, consider the processes Z n defined by the SDE's of the type 
(jO.ip with Z n (0) = and the coefficients a, b replaced by a n ,b n . At last, consider the process Z 
defined by the SDE's of the type (|U.1|) with Z(0) = and the coefficients a, b replaced by a, b. Denote 
f n = Z n (t n ),^ n = zs™ . It is easy to verify that Z n converge weakly in C([0,l],M d ) to Z (see, for 
instance, Proposition 5.1 [H]). Thus, for the sequences f n ,E n , all the conditions of Theorem 13.11 hold 
true with / = Z(t). This means that / possesses a distribution density pf and 

(4.37) sup\p f £(y)-p f (y)\^0. 

y 

Similarly, one can show that, for f n = Z n (t n ), the distribution density pf" exists and 

(4.38) sn V \p r (y)-p f (y)\ - 0. 

v 

Now, ([OSD is provided by the relations (l43TD . (P8j) and 

3*„,tn(yn) = Ps" (Vn ~ %n), Px n ,t n {Vn) = P P {Vn ~ X n ). 

This proves statement (i) of Theorem II. 11 

The proof of (ii') and (iii') can be conducted analogously, with an appropriate changes of the 
truncation procedure and corresponding estimates. Under (B2 exp ), we put, instead of (|4.2p . 

#n(C) = lmax fc < n \\( k \\<SV^- 

By the Chebyshev's inequality, 

P(0 n (C) = 0) < n ■ e-^ 5 ^ 2 < e-' pn . 

with an appropriate p > 0. This and the estimate f|4.34j) provide statement (iii')- 

Under (B2 exp ) and the truncation level given above, the estimates (|4.1(jp - (|4.1|) have their (simpler) 
analogues, and thus the statement of Lemma 14.21 holds. The constant 5 in the definition of the 
truncation level can be made small enough to provide inequality 

d 



n II \ n 



to hold with some no dependent on a, b, U. Then the statement of Lemma 14.31 holds true, also. 
Lemma H3] does not depend on the truncation procedure. Thus the Corollaries 14. 114. 2"1 hold true and 
provide the principal estimate (|4.33p . 
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Under condition (B2 
1 



E 



cxp 



n 



A, a [ X n 



exp;> 

k - 



11 



r=l 



A, b r 



k-1 



n 



3~k-l 



< Cie 



a.s. 



for every A £ K d with some constants C\,C2 dependent on a, b, d, x, Ee^^ 2 . Then the arguments 
analogous to those made in the proof of Lemma 14.21 provide that the estimate (|4.27p holds true for 
every A G M. d . Consequently, the first inequality in Corollary 14.41 holds true for every y > 0. This 
inequality, the estimate (|4.33|) and Theorem 13.11 provide (hi'). This completes the proof of Theorem 
fLU 



4.2. Proof of Theorem ll.2L The implication 2 =^ 1 is obvious. Let us first prove 2 under additional 
supposition (B3). We put, in the notation of Section 3, 6 n {() = 1> f n = ^ Ylk=i 




> cn 



}■ 



n > n* 
otherwise 



with n*, c that will be defined later. Then 

j n d 

Df n = -= ^(Vk)h k =ib r ® e kr , 

Vn fe=lr=l 

where 6 r stands for the r-th coordinate vector in M rf (the proof is straightforward and omitted). 
Since ifi is bounded on U together with all its derivatives, this provides the estimates analogous to 
those given in Lemma [4.21 The statement of Lemma 14.31 is trivial now, since £jj = I-^d (the identity 
matrix in M d ) for every Now, take p = 8(d + l),c G (0, a),?!* > 22+i_ Using Lemma 14.41 we 
obtain the estimate (I4.33|) with the constant C dependent on c, d, tp. The estimate (I4.34p provides 
that P(H n ) < De~ pn for an appropriate D,p > 0. Now the statement 2 follows from Theorem 13.11 

Let us replace the additional supposition (B3) by the condition 1. It is enough to prove the 
statement 2 for n £ mN with some given m £ N. We take m = 2uq with no given in the statement 1 
and have 



d[P„ 



> 



d[p no r , d[p no y 



The function d ^ Pn ° 



d[Pn Q 



d\ d ~ d\ d d\ d ' 
is continuous since A d -almost all points of 



d are a Lebesgue points 

d[P nQ 



In addition, -^^j 



d\ d ls 
R d such 



d\ d d\ d 

(i.e., a points of A rf -almost continuity) for any function g e Li 

not an identical zero due to the statement 1. Thus there exist a > and an open set U C 

Mp lac f/fi 3 ] ac 

that — — * — — > dl(/. Therefore the distribution of + • • • + £ m satisfies (B3). Using what 
we have proved before, we deduce that the statement 2 holds for n G mN, and therefore for n G N. 
This completes the proof of Theorem 11.21 



4.3. Sketch of the proof of Theorem 12.11 We will show that, under conditions of Theorem 12.1 
the following uniform local Doeblin condition holds true. For two measures /ii,U2, denote 

d\i\ dfj,2 



[Mi A/i 2 ](dy) 



mm 



<y), 



("i +V2)(dy). 



Proposition 4.1. For every ball B there exists rig G N, Tb,Jb > such that 



(4.39) 



P X ,T B A P x /, Ts 



> 7b, x,x' e B 
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and, for every n > ng, there exists E ^Z + ,T^ < Tb smc/i i/iat 
(4.40) \p™ T n A P"/ iT »l > IB, x,x' £ B, n> n B . 



Once Proposition 14.11 is proved, one can finish the proof of Theorem 12.11 following the proof of 
Theorem 1 in [4] literally. We omit this part of the discussion and prove Proposition 14. 1| only. 

Proof of Proposition \4-l\ Since (B5) implies (B2 K ) for any k, we can apply Theorem ll.il One can 
easily see that 

i d )> [ wm[£M,<&Jy)]dy. 



pn a pn 

Thus, for any sequence t n — > t > we have, by the statement (i), 



lim inf inf 



pn a p» 



id- 



> inf / mm\p Xjt (y),Px',t(y)]dy. 
x,x'eB j R d 



n-*+oo x ,x'eB I ' b\ 

On the other hand (see [2]), under condition (Bl) the function 

R d x (0, +oo) x R d 3 {x, t, y) ^ p x , t (y) 
is continuous and strictly positive at every point. Therefore, for a given B, 

dcf 1 f 
7s = - inf / mm\p Xtt (y),Px',t(y)]dy > 0, 

Z x,x'£B J^d 

and (|4.39|) . (|4.4U|) hold true for Tg = = 1 and sufficiently large rig. The proposition is proved. 

4.4. Sketch of the proof of Theorem 12.21 It was already mentioned in subsection 2.1 that 
Theorem 12.21 is analogous to Theorem 2.1 [5]. We refer the reader to the paper [5] for the detailed 
proof. Here, we expose a principal estimate only, demonstrating that, in this proof, the truncated 
local limit theorem can be used efficiently instead of the usual one, that was used in [5]. 

Theorem 2.1 [5] is derived from the general theorem on convergence in distribution of a sequence 
of additive functionals of Markov chains, given in the paper [11] (Theorem 1). The characteristics of 
the functionals <p(X),ip n (X n ) are defined by the relations 

f\x) = E[<p°'\X n )\X(0) =x], /*''(*) = E[<f%\X n )\X n (a) = x], s = -, i E Z+, t > s, x E R d . 

n 

The first relation is due to [10], Chapter 6. The second relation was introduced in [TTJ by an analogy 
with the first one. 

The key condition of Theorem 1 [11] is 

(4.41) sup \ft\ x )-f-( x )\^Q, rwoc. 

xeR d ,s=i,te(s,T) 

Here, we need to verify this condition only, since, for all the other conditions, the proof from |5j can 
be used literally. We have 



f^\x) = -F n {x) + - I Fn(y)P: k (dy)=f^- s (x), s<t,xeR m . 

feeN.f <t-s 

We use the decomposition P n = Q n + R n from Theorem 11.11 and write 

f^(x) = -F n (x) + - ]T / F n (y)R^ k (dy) + - £ / F n (y)<Q k (y) dy. 
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The statement (ii) of Theorem 1 1 . 1 1 and the estimates, analogous to the estimates (4.2) - (4.10) from 
[5] , imply that 



sup 

x£R d ,t<T 



fceN,-<t 



F n (y)q: k (y)dy-f\x: 



n 



oo. 



On the other hand, e(n) = | > 1 for k = 6. Thus, by condition (B6) and the statement (hi) of 
Theorem II. 1| 



E 



-F n (x) + - I F n {y)R n k {dy) 



k£N,^<t 

< rT x sup F n (x') 

x> 

uniformly for x G M. d , t < T for any T G 



< n supF n (x') 



1 + E R l 

k<tn 



< 



l + Dn'r . n t + L>^e~ 7fc 

fceN 

This proves (fOTj) . 
References 



n 



oo 



Skorokhod A. V. Asymptotic methods in the theory of stochastic differential equations. - Kiev: Naukova dumka, 
1987. - 328 . 

Win A.M., Kalashnikov A.S., Oleinik O.A. Linear second order parabolic equations. - Uspekhi Mat. Nauk. - 1962. 
- 17, N3. - pp. 3 - 143. 

Konakov V., Mammen E. Local limit theorems for transition densities of Markov chains converging to diffusions. - 
Prob. Theory Rel. Fields. - 2000. - 117. - pp. 551 - 587. 

Klokov S.A., Veretennikov A.Yu. Mixing and convergence rates for a family of Markov processes approximating 
SDEs. - Random Oper. and Stoch. Equations. - 2006. - 14, N2. - pp. 103 - 126. 

Kulik A.M. Difference approximation for local times of multidimensional diffusions. - TViMC. - 2008. - 78, in 



press; preprint is available at |arxiv:ma th. PR/0702175 



Konakov V. Small time asymptotics in local limit theorems for Markov chains converging to diffusions. - 2006. - 
arxiv:math. PR/0602429. 

[7] Ibragimov LA., Linnik V. Yu. Independent and stationary connected variables. - Moscow: Nauka, 1965. - 524p. 
[8] Veretennikov A. Yu. On estimates of mixing rate for stochastic equations. - Teor. Veroyatnost. i Primenen. — 1987. 

- 32. - pp. 299-308. 

[9] Veretennikov A. Yu. On polynomial mixing and rate of convergence for stochastic differential and difference equations. 

- Teor. Veroyatnost. i Primenen. - 1999. - 44. - pp. 312 - 327. 

[10] Dynkm E.B. Markov processes. Moscow: Fizmatgiz, 1963. - 860 p. 

[11] Kartashov Yu.N., Kulik A.M. Invariance principle for additive functionals of Markov chains. - 2006, submitted to 

Teor. Veroyatnost. i Primenen.; preprint is available at arXiv:0704.0508vl. 
[12] Nualart D. Analysis on Wiener space and anticipating stochastic calculus. - Lecture Notes in Mathematics. -1998. 

- 1690. - Ecole d'Ete de Probabilites de Saint-Flour XXV, Springer, Berlin. - pp. 123 - 227. 

[13] Ikeda N., Watanabe S. Stochastic differential equations and diffusion processes. - Amsterdam, North-Holland, 
1981. - 456 p. 

[14] Kurtz T.G., Protter Ph. Weak limit theorems for stochastic integrals and SDE's. - Annals of Probability. - 1991. 

- 19, 3. - pp. 1035 - 1070. 



Institute of Mathematics, Ukrainian National Academy of Sciences, 3, Tereshchenkivska Str., Kyiv 
01601, Ukraine 

E-mail address: kulik@imath.kiev.ua 



